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We use a variational method to construct soliton solutions for systems characterized by opposing 
dispersion and competing nonlinearities at fundamental and second harmonics. We show that both 
ordinary and embedded solitons tend to gain energy when the second harmonic field becomes weaker 
than the first harmonic field. 
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O I. INTRODUCTION 

o ' 

Embedded solitons (ES) represent solitary waves which reside inside the continuous spectrum of a nonlinear - wave 
system. This type of solitons was first reported by Yang et. al. [1] in optical models characterized by opposing 
dispersion and competing nonlinearities at fundamental and second harmonics. More specifically , optical media with 
quadratic x and cubic x^ nonlinear susceptibilities can support ES solutions. The evolution of ES is governed by 
the coupled nonlinear partial differential equations [2] 

^! iu z + ^u 2t + u*v + 7iM 2 w + 472 I v | 2 u = (1) 

and 

iv z - ^Sv 2t +qv+ ^u 2 + 2 72 (| v | 2 +2 I u |> = 0, (2) 

■ where u — u{z,t) and v = v(z,t) represent the fundamental harmonic (FH) and second harmonic (SH) fields respec- 
tively. In writing (I) and (2) we have used optical notations such that z and t stand for the propagation distance and 
reduced time. The quantity —8 is the relative dispersion of SH and 71 2 are the Kerr coefficients. Here q represents 
\D ■ the group velocity mismatch originated by the frequency difference of FH and SH fields. 

Usually, ES's are studied using numerical routines to solve (1) and (2). In view of this one often works within the 
\Q • framework of a simplified physical model where |w| 2 <C \u\ 2 and neglects the cross- phase modulation (XPM) term 
[ (fifth term) in comparison with self-phase modulation (SPM) (fourth term) in (1) . The SPM term in (2) is also 
assumed to be negligible in comparison with its XPM counterpart. Thus we get a truncated model represented by 

iu z + ^u 2 t + u*v + 7iH 2 w = (3) 
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and 



iv z - i(5w 2t + qv + ^u 2 + 4 72 | u \ 2 v = 0. (4) 



For stationary soliton solutions one can use 

u(z, t) = e tkz U(t), v(z, t) = e 2lkz V(t) (5) 

with k, the FH wave number. The partial differential equations of the full model and those of the truncated model 
then reduce to ordinary differential equations given by 

- kU+ -U+ UV + ^U 3 + A l2 V 2 U = 0, (6) 
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- 2kV - ^SV + qV + ^U 2 + 2 72 (y 2 + 2U 2 )V = (7) 

and 

-kU+^U + UV + ^U 3 = 0, (8) 

-2kV -^SV + qV + -U 2 + A l2 U 2 V = 0. (9) 

Here the dots denote differentiation with respect to t . Linearization of the equations in (1) and (2) [Full model] as 
well as in (3) and (4) [Truncated model] shows that both models support ordinary soliton solutions in the regions 

0</c<|if<5>0, k> max |o, || if S < (10) 

and embedded soliton soliton solutions in the regions 

k > max jo, !} if 5 > 0, < k < | if S > 0. (11) 

The object of the present work is to derive a straightforward analytical model for comparing the properties of 
soliton solutions supported by the pair of equations representing the full and truncated models .In doing so we shall 
consider the cases of ordinary and embedded solitons separately . To achieve this we shall envisage a variational 
approach to the problem, where one begins with a Lagrangian for the system under consideration and constructs the 
so-called effective Lagrangian by taking recourse to the use of trial functions for the field variables. Understandably 
, the trial functions will involve a number of unknown parameters . As we shall see the effective Lagrangian will 
provide a natural basis to determine these parameters . 

In the above context we note that (6) and (7) , resulting from the full model , follow from an action principle 
. In contrast to this , (8) and (9) pertaining to the truncated model are non-Lagrangian . But the latter set of 
equations are based on physically founded assumptions . This led Kaup and Malomed [2] to adapt the variational 
approach to the seemingly flawed system represented by (8) and (9) . In their method one starts with the Lagrangian 
of the full system and drops the term containing V to construct an expression for the effective Lagrangian by using 
the trial functions for U and V . Further , the implementation of the Ritz optimization procedure to evaluate the 
variational parameters requires one more approximation . Wc claim that the results in Ref (2) can be rederived and 
reexamined without taking recourse to the use of this two - tier approximation .In particular , we find that if we work 
with the effective Lagrangian of the full system , construct equations for the variational parameters and then use the 
approximation U <C £7 ,we automatically arrive at the results of Kaup and Malomed . More significantly , the method 
followed by us provides a natural basis to examine how the results for U and V for the full model differ from those 
of the truncated model . One of our main objectives in this work is to compare the results of the full and truncated 
models and thereby gain some physical weight for the problem . 

We begin section II with the Lagrangian of the full system and construct the expression for the effective Lagrangian 
using some trial functions for U and V. We then apply the Ritz optimization procedure to obtain equations for the 
parameters of the trial functions and examine how the results of Rcf. 2 arc obtained for V <C U . In section III we 
compare the results of U and V for the full model with those for truncated model. We represent the results for both 
ordinary and embedded solitons. 

II. VARIATIONAL FORMULATION OF (6), (7), (8) AND (9) 

A. Lagrangian representation 

Our analysis for the properties of ordinary and embedded soliton solutions supported by the full and truncated 
models will involve essentially a Ritz optimization procedure [3] based on the variational functional for (6) and (7) 
. It is easily seen that these initial-boundary value problems can be converted to a variational problem with the 
Lagrangian written as 

L = J [ - kU 2 - (2k - q)V 2 - l -U 2 + 5 -V 2 + U 2 V + 

^±-U 4 + A l2 U 2 V 2 + l2 vAdt. (12) 



FIG. 1: U and V as a function of t for non - embedded solitons. 



In the Ritz optimization procedure, the first variation of the variational functional is made to vanish within a set of 
suitable chosen trial functions. We thus introduce the ansatz [2] 

U = Asech(V2kt) and V = Bsech 2 (V2kt) (13) 

for the time - dependent parts of the FH and SH fields. Here the amplitudes A and B are variational parameters. The 
inverse width \[2k will, however, not be varied. Inserting (13) in (12) and carrying out the time integral we obtain 

(L) = f - 4kA 2 - 2(2k - q)B 2 + -SkB 2 + 2A 2 B + 

7l A 4 + | 72 A 2 S 2 + | 72 B 4 ], (14) 

the effective Lagrangian for U and V in (13) . The Lagrangian in (14) represents a specific function of the parameters 
only. Optimization with respect to parameters will yield a system of equations which when solved will determine U 
and V within the chosen set of trial functions and a concomitant approximation for the true solutions. This is the 
route we follow to determine the values of the parameters A and B. 



B. Variational parameters and truncated model 

From the vanishing conditions of and ^j^- we obtain 

-2k + B + ^A 2 + ^- l2 B 2 =0 (15) 
5 

and 

OO Q Of; 

A 2 + —-i2A 2 B - 2{2k - q)B + UkB + ^ 72 B 3 = 0. (16) 
5 5 35 

Understandably, these equations determine the parameters of the full model. To go over to the truncated model we 
can choose B<i and neglect B 2 and B 3 in (15) and (16) to get 

-2fc + B + 7 iA 2 = (17) 

and 

32 8 
A 2 + —-f 2 A 2 B — 2{2k - q)B + -SkB = 0. (18) 
5 5 

These equations were obtained by Kaup and Malomed [2] first by neglecting the last term in (12) and then again 
neglecting the contribution of the term 47 2 C/ 2 I /2 while taking variation with respect to A. But we have shown that 
this type of two - tier approximation is not essential to make a transition from the full to the truncated model. 
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FIG. 2: Z7 and V as a function of t for embedded solitons. 



III. SOLITON SOLUTIONS 

We have seen that when the wave number k falls into the region (10) both full and truncated models have ordinary 
soliton solutions . To see how , in this case ,the results of U and V for the full model differ from those of the truncated 
model we have chosen to work with k = 0.25, 71 = —0.05,72 = —0.025,(5 = l,q = 1 . As for the full model , we 
use these values in (15) and (16) to get three values for B , namely B 1 = 0.4824, B 2 = 5.6895, B 3 = 19.4379 . We 
find that A values corresponding to B x and B 3 arc imaginary while A becomes a real number equal to 7.2109(Af) 
when calculated by using the value of B^iB^ ) .The corresponding results for the truncated model are B^ = 6.7227 
and A \ = 11.1559 . The superscripts F and T refer to the full and truncated models . We shall also use similar 
superscripts on U and V .The ordinary or non-embedded soliton solutions arc shown in Fig. 1 . From this figure it is 
clear that the curves for U T and V T are more peaked compare to the curves for U F and V F . It will , therefore be 
interesting to examine how the behaviour of U F , V F , U T and V T is affected in the case of embedded solitons . 

In consistent with (11) we take k = 0.6963, 71 = —0.05,72 = —0.025,(5 = l,q = 1 for the embedded soliton . In 
this case we find B F = 4.9450, A F = 5.6500, Bj = 6.3822 and A\ = 9.9896 . In Fig. 2 we display the curves for 
embedded solitons . In this case also the curves for U T and V T are more peaked than the curves for U F and V F . 
But looking closely into the curves in Figs. 1 and 2 we see that in the case of embedded solitons the curves for U T 
and V T fall off more rapidly than their non-embedded counterparts . 

In the model considered in this work the energy of the soliton is given by E = f*°° (|w| 2 + 2|w| 2 ) dt. From the 
results in Figs. 1 and 2 it is clear that E T > E F for both ordinary and embedded solitons. Understandably, E T and 
E F stand for the soliton energies obtained by using the truncated and full models. It will, therefore, be an interesting 
curiosity to verify how the approximation V<(7 affects a typical experiment. 
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